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■ We develop a gauge-invariant linear response theory for relativistic Bardeen- 

^ ■ Cooper-Schrieffer (BCS) superfluids based on a consistent-fluctuation-of-the order- 

parameter (CFOP) approach. The response functions from the CFOP approach 
satisfy important generalized Ward identities. The gauge invariance of the CFOP 
\ theory is a consequence of treating the gauge transformation and the fluctuations 

of the order parameter on equal footing so collective-mode effects are properly in- 
1^^- eluded. We demonstrate that the pole of the response functions is associated with 

the massless Goldstone boson. Important physical quantities such as the compress- 
fH ! ibility and superfluid density of relativistic BCS superfluids can also be inferred from 

\ our approach. We argue that the contribution from the massless Goldstone boson is 

crucial in obtaining a consistent expression for the compressibility. 
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I. INTRODUCTION 



The Bardeen-Cooper-Schrieffer theory of fermionic superfluids (see [1] for a review on its 
apphcation to conventional superconductors) has provided not only an explanation for con- 
ventional superconductors and other fermionic superfluids but also a paradigm for studying 
macroscopic quantum coherence due to interactions. Soon after its discovery, the challenge 
of how to cast its formalism in a gauge invariant form when a charged system is interacting 
with an electromagnetic (EM) field has drawn broad interest. At the linear- response-theory 
level, there have been two major approaches. Nambu in his seminal paper [2] formulated 
this challenge in terms of generalized Ward identities (GWIs) and proposed an approach 
based on a set of integral equations for finding a gauge-invariant EM vertex which governs 
the kernels of response functions. This approach has been reviewed in Schrieffer's book on 
superconductivity [H] and also finds applications in other branches of physics such as nuclear 
matter interacting with neutrinos j^. 

There is another approach based on consistent fluctuations of the order parameter 
(CFOP), which is the main theme of this paper. In this approach, the effects of gauge trans- 
formation from the gauge field are balanced by the fluctuations of the order parameter in a 
consistent fashion. This is made possible by treating the terms induced by the gauge field 
as well as the fluctuations of the order parameter equally in the perturbative Hamiltonian. 
Although the kernels of response functions from this approach are not a solution of Nambu's 
integral equations, the CFOP formalism is manifestly gauge invariant and we will show that 
GWIs are satisfied. This is because being a solution of Nambu's integral equations is only 
a sufficient condition for satisfying the GWIs, but not a necessary condition. Importantly, 
this approach reproduces the compressibility correctly as that from the equations of state 
and this demonstrates self-consistency of the CFOP approach. 

The theory of CFOP has an interesting history. Kadanoff and Martin [3] first proposed 
this approach in a less complete form by considering only the phase fiuctuations of the 
order parameter. Their idea is to decompose the three-particle Green's function in a way 
that can respect gauge invariance. Betbeder Matibet and Nozieres [HI and Kulik et al. 
[gI independently formulated this approach in more complete forms with both phase and 
amplitude fiuctuations later on. This approach has also been formulated by the Keldysh 
formalism with time-ordered Green's functions in Ref . 0] • After its successful application to 
conventional superconductors, a generalization of this approach to nonrelativistic fermionic 
superfiuids satisfies important sum rules and has been applied to ultra-cold atomic systems 



10|. It has also been discussed in the physics of neutron stars pJLj]. Here we base on our 
formalism of a relativistic version of the BCS theory [l2| with Kulik's approach to CFOP and 
develop a gauge-invariant linear response theory of relativistic fermionic BCS superfiuids. To 
demonstrate the versatility of this approach, we will address the collective mode associated 
with the massless Goldstone boson in the symmetry-broken phase, the density susceptibility 
which governs the compressibility, and the superfiuid density. 

To further contrast these two approaches, we also present the relativistic version of 
Nambu's integral equations for the EM vertex. There have been attempts to find an it- 
erated solution based on the random phase approximation (RPA) [sj. In nonrelativistic 
BCS superfiuids it is possible to argue that the RPA-based theory satisfies the correspond- 
ing GWI. To our knowledge, neither a relativistic version of Nambu's integral equations nor 
a complete proof of the gauge invariance of the RPA-based linear response theory for rela- 
tivistic BCS superfiuids have been explicitly presented. Since a major goal of this paper is 
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to advocate the CFOP theory of relativistic Fermi superfluids, we will limit our discussions 
on Nambu's integral-equation approach. 

This paper is organized as the following. Sec. [TTl briefly reviews a microscopic theory for 
relativistic BCS superfluids that will be the foundation of this work. Sec. IIIII presents the 
CFOP formalism and we explain in more detail how our theory respects gauge invariance 
in Sec. IIVI Sec. |V] gives the explicit expressions of the response functions from our CFOP 
approach. Sec. IVll Sec. I VI II and Sec. IVIIII show some applications of the CFOP theory to 
the collective modes, compressibility, and the Meissner effect for a relativistic BCS super- 
fluid. We briefly discuss a relativistic version of Nambu's integral-equation approach and its 
associated GWIs in Sec. [IXl Sec. |X] concludes our work. Some details and conventions are 
given in the Appendix. 



II. MICROSCOPIC THEORY OF RELATIVISTIC FERMI SUPERFLUIDS 

Several relativistic models of a two-component BCS superfluids have been formulated in 



Ref s. [I2l 4l7l] and we briefly review the model following the BCS-Leggett mean field theory 



I2I . Il5| without any external gauge field here. The Lagrangian density is 



£(x) = J2 M^l^d^ -m + /i7°)V^. + Cii^), (1) 

where ip, ip are Dirac spinors which denote the fermion fields with mass m and chemical 
potential /i. The representation of the 7— matrix and some useful properties are given in 
Appendix \M Throughout this paper, we take the convention e = c = h = 1 and use a to 
denote the pseudo-spin '\ and I with "1^= — I and a = —a. The pseudo-spin may refer to 
some internal degrees of freedom such as the color indices in quantum chromodynamics. £/ 
describes the attractive pairing interactions between particles with different pseudo-spins 



and it takes the form 13 



£,(x) = ^?(^^^«75^J(^4^75C^^t ), (2) 

where g is the attractive coupling constant, and the charge conjugation matrix C is defined 
as C = i'jol2- The gap function that is also the order parameter is given by A(x) = 
gi^p'^Ci'^^ip]). The standard BCS approximation then gives 

C-Bcsi?^) = Yl Mn^d, - m + /i70)V'^ + A*(V'?^«75^J + A{^^zj,Cij^). (3) 

The corresponding form of the Hamiltonian density is then given by 

Hbc5(x) = ^-(-^^ ■ V + m - /i7°)V'<. - A*(^^^«75^4) - A(^^Z75C^^^). (4) 

In the broken-symmetry phase the order parameter may be chosen to be real. Here we 
present our theory in Matsubara formalism, which is applicable to both zero and finite 
temperature T. We will focus on T = results and a generalization to finite T within the 
BCS approximation is straightforward. To simplify the notation, we group the imaginary 
time T = it and x as a four- vector x = (r, x) and define 

0{x) = e^^^^"C(x)e-^^^sr (5) 
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where Hbcs = J d^^T^BCsi'^)- The single particle Green's function and anomalous Green's 
function are given by 

G{x,x') = -{TrMx)i^^{x')]), F{x,x') = -{TAM^)^Ji^')C]), (6) 

where denotes the r-order of operators. When Hbcs is time- independent, G and F 
depend only on the difference r — r'. Let x+ = (r + 0+,x). The gap function can be 
expressed as 

A(x) = (7Tr[275F(x,x+)]. (7) 

The number density for each species is defined by ^^-(x) = {ijj„{-x)'y'^ipa{x.)) = {ipl{x.)ilj„{x.)). 
It can be also calculated from the single particle Green's function no-(x) = Ti^j^G^XjX 
Therefore, the total fermion number is given by 



n 



J £x{n-^{x) + ni{x)) = 2 J £xTi [7°G(x, x^)] . (8) 



It can be shown that the Green's function and anomalous Green's function satisfy the 
following equations of motion 

[-7° 9, + . V - (m - /i7°)]G'(x, x') + i^i^F{x, x') = 5{x - x') l4x4, (9) 
[-7°^^ + . V - (m - /i7°)]F(a;, x') + iM^G{x, x) = 0, (10) 

where G{x,x') = GG^{x',x)C, F{x,x') = 7°F"''(x', x)7°, l4x4 is the four- dimensional iden- 
tity matrix, and 6{x — x') = 6{t — r')(5(x — x'). Here we define the fermion four-momentum 
at finite temperature as P = {ioJn, p), where a;„ is the fermionic Matsubara frequency given 
by iJn = (2n + l)TTkBT, where fc^ is the Boltzmann constant. The quasi-particle energies 
are given by = with = ep ± /i and ep = a/p^ + m^. With the help of the 

energy projectors 

A,(p)4|l± -'°'^"P + '"' ]. (11) 
the solution of G and F in momentum space are 

F(P,/i) = [-^ — : — ^— + — : — «75, 13 

lUJn - E„ lUJn + E^ lUn - E+ lUJn + E+ 



where = -^(1 + -^) andf^^ = — The energy projectors A+(p) and A_(p) project 

out the contributions from the fermion and anti-fermion respectively. In the nonrelativistic 
limit where |p| ^ m, |/x — m| <C m and A ^ m, one gets A+(p) ~ 1 and A_(p) ~ 
so our expressions reduce to the well-known nonrelativistic results. By taking the Fourier 
transform of Eq.([H]), the fermion number is given by n = 2 Tr[7°G'(P, /i)], which is the 
number difference between the fermions and anti-fermions 



n = n+ — 
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where n± denote the density of fermion and anti-fermion. The Fermi momentum kp is 
defined by n = 2kp/ (Svr^), and the Fermi energy is ei;' = ^ykp + m?. The Fourier transform 
of Eq. ([7]) gives A = (7 Tr[i75F(P, /i)], which leads to the gap equation 



p 



E- 



+ 



(15) 



In the nonrelativistic limit, the number equation ( fT4|) and gap equation ( ITSl) reduce to the 
well-known nonrelativistic results except a factor of 2 on both right-hand sides. This factor 
of 2 comes from the fact that we have introduced the pseudo-spin, which brings two times 
more degrees of freedom. 

Since the model is not renormalizable in 3 + 1 dimensions, a regularization or a momen- 
tum cutoff A is needed. The relativistic limit of the BCS state depends on the Compton 
wavelength Ac = Tn~^ 15|. li kp 1 / \c = the system evolves into the relativistic regime. 



It has been shown that this model can be generalized to describe the BCS - Bose-Einstein 
condensation (BEC)- relativistic BEC crossover 15| of Fermi gases. 

To further compactify our expressions, we reformulate the relativistic BCS theory in the 
Nambu formalism pL,, |2i]. This is more convenient for the discussions on the linear response 
to an external electromagnetic field. We introduce the Nambu-Gorkov spinors 



(^t(x), ^r(x)c). 



Moreover, we define 



-((Ti +^a2),(T- 



- 0-1 
2^ 



cr3),cr- = ^(ctq - CTs)- (16) 



in Nambu space. One can show that A = gi^^i'-y^a-'^) and its Fourier transform is Aq 
(7 ^p(^pZ75(T_^p+q). Similarly one can show that 



C-Bcs = ^{il^d^ - m + /x7o(T3)^ + ^(Az75(t+ + /\*i-i^aJ)m 



(17) 



The Lagrangian density can be written as Cbcs = where the inverse propagator in 

momentum space is given by 



G~\P,i2) = -(T,[^p^p]) = (iw„ + /ia3)7°-7-p-m + Az75ai. 



(18) 



After evaluating the inverse of the right-hand-side, one gets the expression of the propa- 
gator (see Appendix [B]) 



+ 
+ 



iur, — E~ 



+ 



2 M+2 



u. 



+2 



,+2 



+ 



iUn — E^ 
A+(p)A 

{tUnf - E-^ 



+ 



-)A+(p) + ( 



u 



-2 



A-(P)A 

[zUn)' - E+^ 



+ 



+ 



iu„ — E'^ 



-)A-(P) 



,-2 



iUn + E^ iUn - E. 



-)A-(P) 



7°a+ 



^750-+ + 



A+(p)A 



+ 



A-(P)A 

[^U^n? - E-^ 



i75Cr_. (19) 
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From Eqs. ( !T2|) and ( !T3|) . one finds that 

^^(P^_^) G(P,-/i) ) 



(20) 



Moreover, the number equation (fT4l) and gap equation ([71) can also be rewritten in Nambu 
space as 



n 



(21) 



The expression f ll9p of the propagator in Nambu space can be further simphfied to a more 
instructive form. We define the operator Ep = 7''(7 ■ p + fn) — ficr^ — A'y^i'y^ai in Nambu 
space and also introduce the projectors 



A+(P) 



A+(P) 
A.(p) 



A_(p)^l-A+(p) 



A-(P) 
A+(p) 



Then the propagator becomes (see Appendix lUj) 



^iun - -Ep iUn + E~ iUn + E+ lU^ - E+ 



where the coefficients are given by 



(E^Ti^p)A^(p) 
S = ^ '^P 



(E±±4)A^(p) 



2E± 



(22) 



(23) 



(24) 



Those coefficients are the counterparts (in Nambu space) of the coefficients and v^. 



III. GAUGE-INVARIANT LINEAR RESPONSE THEORY 

We consider fermions and anti-fermions coupled to a weak external EM field A^(x). The 
derivative 9^ in the Lagrangian density ([T]) should be replaced by the covariant derivative 
= + iy4^(x), which results in the interacting term £a(x) = — X]o-=t4, '^o-^^'^o-^m ~ 
—J^A^. In Nambu space, one can show that Ca = —^^'^cr^.'^A^. 

The corresponding Hamiltonian density is 'Hyi(x) = — £^(x). Gauge invariance of a 
microscopic linear response theory with respect to an external EM field is made possible by 
considering the perturbations due to the fluctuations of the order parameter in a consistent 
fashion. The nonrelativistic version of the CFOP method has been extensively studied 
and here we will develop a relativistic version of this method. 



In equilibrium, the order parameter is given by A. We assume that the deviation of the 
order parameter from its equilibrium is small and denote the small perturbation by A'(a;). 
Therefore, A in Eq. fll7l) is replaced by A — > A + A'. Then the Hamiltonian density splits 
into two parts: the equilibrium expression and the part containing the deviation. Explicitly, 
T-Lbcs = T-Lbcso + T-i' where 

Ubcso = ^(-i7-V + m- /X7V3 - Ai75ai) ^ (25) 
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and 

n' = ^ ( Aii75(Ti + A2Z75(^2 + 4(^3} ^. (26) 

Here A' = — (Ai — iA2) and A'* = — (Ai + iA2). Ai and A2 are the negative real and 
imaginary parts of the fluctuations of the order parameter. The Hamiltonian Hbcs becomes 

Hbcs = Hbcso + H' (27) 

= ^ ^^p(7 ■ P + - Ai7°0-3 - Ai75Cri)^p + ^ l|'p+q(AiqZ750-i + A2qi75Cr2 + ^qCTsj^p 

p pq 

The interaction term may be considered as a scalar product $q ■ S, where 

#q = (Aiq, A2q, A^q)^, S = ((TiZ75, ^2^75, ^^37'')^, (28) 

are the generalized external potential and the generalized vertex function. To calculate the 
linear response of a relativistic Fermi superfiuid to the perturbation H', we introduce the 
response-function vector fj: 

r/(r,q) = ^(^p(r)SM/p+q(r)), (29) 
p 

where 773 = corresponds to the current due to the external field and 771,2 denote the 
perturbations due to the fluctuations of the gap function. The covariant index should not 
be confused with the chemical potential. The linear response theory is then written in a 
matrix form 

r/(r, q) = Q{t, q) ■ $q 

Qn{T,q) Qi2(r,q) Qi3{T,q) 

Q2i(r,q) g22(r,q) Q^r^cO 1 I A2q I . (30) 





r Alq\ 




A2q 







The response functions Qij are 

g,,(r-r',q) = - ^(T,['l'p(r)S,xI/p+q(r)^p,+q(r')S,^p.(r')]). (31) 
pp' 

Using a Fourier transform and making use of Wick's theorem, we obtain 

Qij{i^li, q) = TrT^ ^ ^iGp+^^p'+^iiUn + i^i)tjGp^p>{iuJn) 

iuJn pp' 

= TttJ2Y1 ^^^(^ + ^' ^')^AP^ (32) 

iuln P 



where Q = {ifli,q), Qi is the boson Matsubara frequency, Gp^p'^iUn) = Gp{i'^n)Sp,p', and 
Gpitun) ^ G{P,ix) = -1 
response matrix is given by 



GpUun) = G{P,ix) = - — ^r7°- Inserting the above relations into Eq. (15^ . the hnear 

ILiJn — -C/p 



Q,,{^Q,,cD=TJ2J2^4^^- — ^7°). (33) 

p ^ lUJn+ till - -Ep+q lUn - Ep ^ 
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Next we show that if the fluctuations of the order parameter are formulated as shown in 
Eq. (126|) . our microscopic hnear response theory is exphcitly gauge invariant. Applying the 
condition 771 2 = — |Ai^2 consistent with the gap equation to Eqs. (l30l) . we flnd 

QIIQ22 — Q12Q21 

Q11Q22 — Q12Q2I 

where Qu = f + Q22 = g + ^22- After substituting the results into 

= = g^i Ai + g^2A2 + g^3 (35) 

we get = K^"'A^ = {K^" + 6K^"')A^. Here fsT^^ = Q^g and 

SK^" = ^11^32<?23 + Q22Q31Q13 - Q12Q31Q23 ~ Q21Q32Q13 ^3g^ 

Q11Q22 ~ Q12Q2I 

The gauge invariance condition q^J^ = leads to q^K'^'^^Q) = 0, where q^ = Q = {iQi,q) 
is the covariant form of the four-momentum. This condition is explicitly satisfled if the 
response functions satisfy the generalized Ward identities 

Qf^Qsi = -2'iA(52i, 

QtJ.Q32 ~ — 2'iA(522; 

g^Q^3^ = -2zAg^3. (37) 

The derivations of these GWIs will be given in a moment. Firstly we show that the response 
functions indeed satisfy those GWIs. We observe that 



Q11Q22Q23 + Q22Q21Q1S — Q12Q21Q23 " Q21Q22Q 



V 

13 



Q11Q22 — Q12Q2I 

= -22AQ^3 + 2iAg^3 = 0. (38) 

The proof of the GWIs (1371) is sketched here. Our starting point is the expression fl32l) . 
In what follows we will use the covariant form of the four- momentum = P = {iun, p) and 
P interchangeably. Moreover, = TJ2iu;„ J2p- If '^^ apply the analytical continuation 
icOn CO + i5, then = J By using the notation p = 'j^icon — 7 ■ p we can express 

the bare and full propagators in Nambu space as 

Go^iP^ ^j) =f-m + cr3/i7°, 

G-\P,ii) = G^\P,li)-% (39) 

where the self-energy in Nambu space is S = — Ao"iZ75 and should not be confused with the 
effective vertex function S. Those expressions give 

a3G-\P + Q,ii)- G-\P, fi)a3 = fa^ + 2zAa2Z75. (40) 

Eq. (UnD will lead to the GWI ([ST]). One can show that 

g^g^i + 2zAg2i 
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= Tr J] [{a,iG-\P + Q,fi)- G-\P, fx)a,) G{P + Q, ii)a,t^,G{P, /i)] 
p 

= -2Tr5^[a275G(P,/i)] 
p 

= 0, (41) 
where Eqs. ( lElj) has been used. Similarly, for the second GWI we can show that 

22 

= Tr J] [a3a2^75G(P, /i)] - Tr J] [G{P + Q, (1)115^2^3] 
p p 

= -22Tr5^[aiZ75G'(P,/i) 



4? 

= --A, (42) 
9 

where in the last line we have used Eq. (pTI) . Therefore we get the second GWI qfiQ32 — 
-2iA{Q22 + I A) = -2iAQ22. For the last GWI of Eq. §7^, we have 



g^Qg + 2zAQ^3 

= Tr ^ [YG{P, fi)] -TtJ2 [GiP + Q, f^)r] 

p p 
= 0. (43) 

This completes the proof. 

IV. MORE ABOUT GAUGE INVARIANCE 

We have seen how the gauge invariance condition is satisfied by our response functions. 
Here we will clarify some subtleties from a generalized interaction picture. In our linear 
response theory, the Lagrangian density after the BCS approximation is given by 

^BCS = ^BCSO + ^' 

= ^ {i-f^df, -m + fi-f'^a-s + Aij^ai)"^ - ^ {Aiij^ai + A2'i75f^2 + ^(^s) ^- (44) 

One can show that it is invariant under the generalized infinitesimal gauge transformation 

^ ^ (l+i(J3X)*, 

A ^ A, 
Ai ^ Ai, 

A2 ^ A2 + 2Ax. (45) 
Under this transformation the two parts of the Lagrangian density transform as 
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= C-BCSQ - ^(tJx^ + 2x^Ai75a2^, 
a ^ C - 2x^'Ai75a2^ + ^a^$x^- (46) 

Therefore Lbcs is invariant under the generahzed infinitesimal gauge transformation fH5|) . 
Here we emphasize that the transformation f H5|) only keeps terms linear in x in the linear 
response theory. 

The original mean-field Lagrangian density ( [T7|) with a real A is not gauge invariant if 
the order parameter A is not perturbed by the gauge transformation. To circumvent this we 
assume that the effects are absorbed into the fluctuations of the order parameter A' while 
the equilibrium value of A is unchanged. From A(x) = g {ip'j^ Ci'^^ipi) ^ind A2 = ImA' we see 
that only the imaginary part of the order parameter is perturbed by ip^i — )■ {l + ix)ipu while 
the (negative) real part is not, Therefore the perturbations from the gauge transformation 
are 5Ai = and 6A2 = 2ixA. It is the gauge transformation of A2 that cancels the term 
associated with A in the generalized gauge transformation and leads to the gauge invariance 
of the Lagrangian density. 

As shown in Section JITTI the perturbative Lagrangian density can be written as ^^^^ ■ I]\E', 

where $ and S are deflned by Eq. f l25]) . In fact ^ may be viewed as the generalized external 
gauge fleld and S may be viewed as the generalized vertex function. In our theory, there 
are three different spaces: (i) the two-dimensional Nambu space where the Pauli matri- 
ces live, (ii) the four-dimensional representation space of the Clifford Algebra in which the 
7— matrices live, and (iii) a three-dimensional space which we will deflne as the generalized 
gauge space, where the generalized external potential $ and generalized vertex function S 
are deflned. Thus the transformation (145|1 is the corresponding generalized gauge transfor- 
mation of the generalized external gauge fleld. Explicitly, 

. . / \ 

* ^ $ + 2Ax . (47) 

V -9,x J 

In momentum space —d^x becomes —iq^X- We deflne the generalized external momentum as 
q = (0, 2iA, q^Y^ in the generalized gauge space. Then the generalized gauge transformation 
f H7|) can be written as 

iqx- (48) 

We saw that the Lagrangian density is invariant under the generalized gauge transformation 
f H5|) and now we want to flnd the corresponding GWI. By using the form of generalized 
external momentum, Eq. fHOj) becomes 

a^G-\P + Q,fi)- G-\P, fi)a3 = + 2zA(T2Z75 = ■ S. (49) 

This is the GWI associated with the generalized gauge symmetry. 

Next we adress what the conserved generalized current associated with this gauge trans- 
formation should be. Using the self-consistent condition 1712 = — ^Ai2, Eq. fISU]) can be 
written as 

\ /Qn Q12 Q'is\ /Ai\ 

] = { Q21 Q22 Q'23 \ \ ]■ (50) 
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We then define the generahzed current J = (0, 0, J^Y' and three generahzed response- 
function vectors 

Ql = I I , Q2 = I ^22 I , = I I • (51) 





Q33 

Then the generalized current (150|) becomes 

J = (Qi,Q2,Q^)-^, (52) 

The GWIs fl37|) for the response functions can also be written as 

q^-Qi = fori = 1,2, 3. (53) 
Thus the GWIs directly lead to the conservation of the generalized current 

■ J = (q^ ■ Ql, q^ ■ Q2, q^ • Q^) • $ = 0. (54) 

This gives 

g^J'^ = 0. (55) 

Therefore J is indeed the conserved current associated with the generalized gauge transfor- 
mation. We see that the generalized gauge transformation leads to the usual U{1) gauge 
invariance of our linear response theory. Importantly, the GWI ( H9|) for the generalized 
vertex function is exact since there are no high order corrections to the vertex in the linear 
response theory. Thus we have proved that our CFOP theory is indeed gauge invariant. 



V. EXPRESSIONS OF THE RESPONSE FUNCTIONS 

It will greatly simplify our expressions of the response functions from the CFOP approach 
if we sort out the odevities of them first. Here we list the main results and leave the details 
in AppdexjEl The odevity of the response functions about the four-momentum Q = {ifli, q) 
is 

QijitQi, q) = {-lf''+''')Q,^{-^Ql, -q). (56) 

The odevity about the gauge indices i,j is given by 

Qjiiin^, q) = (-l)(^"+^'^)g,,(^^],, q). (57) 

The odevity about the boson Matsubara frequency iQi is relatively complicated. For i = 
j = 1,2 we have 

Qii{iVli, q) = Qiii-i^i, q). (58) 

For i = j = 3 we have 

(^^^^ q) = ( ^^E;'^^ ^\ ^ Y ^ ° °' ^ ^ ' ^ 'n (59) 
^•^•^^ [_ —Q'^^{—iili,q) it fi = 0, u = I or fi = I, u = ^ ' 
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For i = 1, j = 2 we have 

gi2(in;,q) = -gi2(-A,q). (eo) 

For i = 1, j = 3 we have 

«w^o,.<,)^{_f(-f^;;'^)-f;::° (61) 

For i = 2, J = 3 we have 

After sorting out the odevities, the expressions of the response functions can be derived. 
After summing the Matsubara frequencies, Eq. ( l33l) becomes 

Q,jitni,q) = Y, J de, J de2Tr (^ -^^^'^ ~ [^^^ t^S (e, - E^+^)^%5{t2 - 4)7°)- (63) 

The 5-function operator can be decomposed as (see Appendix ICj) 

5{t - Ep) = Wp5(e -E^)+ v^5{t + E^) + upie + E^) + t)+5(e - E;). (64) 
We define the coherence coefficients as 

(u^v^),, = Tr[E,w;+^7°S,^)j7% (u^v^h = 

(v^v^),, = Tr[E,t)±+q7°S,^)j7°], (v^v^h = Tt[±,v;^^^%v^^% (65) 
The response functions can be exphcitly written down as 



Qij{iQi,q) 



+ 



{f{E-^^)-f{E-)){u-u-: 




(1 


- /(^p+qy 


^-f{E-)){u- 


~v 


)ij 


Ep+q - Ep - i^i 








+ Ep - tQi 






(l-/(Ep^J-/(E+))(n 




%^ 


(/(^p+q: 






')■■ 
Jtj 


+ E+- zQi 






Ep, 


^^-E+- zQi 






{l-f{Ep^^)-f{E^)){v- 






(/(^P+q) 


'~f{Ep)){v- 




)ij 


^p+q + ^p + ^^l 






Ep^ 














- /(^p+q) 


'~fiE;))iv~ 


v+ 


)ij 








-E'p+q 


+ E+ + iQi 






{l-f{E;^^)-fiE-)){u 






(/(^p\q: 


^~f{E-)){u^v- 


)ij 


^p+q + ^p 








-q - ^p + T'^l 






{fiE;^^)-fiE;))iuv: 




(1 


- /(-£^p+qy 


^~f{E;)){u^ 




)ij 


^p+q -E+ + t^l 






^P+q 


+ E+ + tVli 
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+ 



ifiK^^) - f{E-)){v+u-\, (1 - f{E;^^) - f{E-))iv+v-) 



^P+q + ^P - 



Ep - in 



E;^^ + E+-ini 

If Qij is an even function of iVLi, the expression reduces to 



{m^^)-m)){E^^^-E-) 
(^p+q-^p)'-(^w 

(1 - /(^p+q) - fjE^)) {E-+C, + ^p-) 

(^p-+q+^p)'-(^w 

(1 - /(i^p+q) - /(^p+)) (^p+q + 



(/(i^p+q) - m)) {E^^^ - e;) 

+ terms with super-indices (+ ^ — 



+ 



If Qij is an odd function of iQi, the expression reduces to 



+ 



/(^P+q) - fjE^) 

l-/(^P+q)-/(^p: 
{E^^^ + E^f - [tni] 

{E^-.^-Eif-mi)- 

terms with super-indices (+ 4-)- 



■{{u V -{v u )y) 



■(( 



U U )ij — [V V 



(66) 



(67) 



(68) 



The coherence coefficients are listed in Appendix jFl 



VI. COLLECTIVE MODE CONTRIBUTION TO RESPONSE FUNCTIONS 

In Section IIIII one saw that the conservation of the induced current J'^ = requires 
q^K'^'^{Q)A'^ = 0, where A'^ and K'^^[Q) are an external field and its corresponding response 
kernel. In general, if one applies a gauge transformation = A^ — iq^Xi the response kernel 
transforms as K''^^{Q) K^^^^Q) and the gauge invariant condition requires q^K^^''{Q)Ai, = 
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0. From the gauge invariance condition q^K'^^ {Q) = 0, one obtains ix 
Using the gauge invariance condition again one obtains 



qa 



113 



(69) 



This is the general expression of how the response kernel transforms under gauge trans- 
formations. Clearly, the zeros of q^^K'^'^q^ = 0, if exist, indicate the presence of collective 
excitations. In the following we will show that the massless Nambu-Goldstone boson indeed 
contributes. We simplify the expression of K^'^ by the response functions we obtained so 
far. Following the discussions in Ref. 



we rewrite K^'^ in a more compact form 



with 



-"■0 



^533 



13 



Q 



11 



Q32Q'23 
Q'2 



22 



where 



V32 ~ V32 ^ V31' V22 ~ V22 ~ 

Vll Wll 
By using the GWI (1371) and the symmetry (]E9p . one can show that 



(70) 



(71) 



(72) 



K\xqsKt 



lap 



qp 



(73) 



This implies that K^'^ can be thought of as a functional of Kq^ and the denominator of the 
second term is expressed as q^K'^^q^ = A^t^Q^i- This expression is useful in the analysis 
of the dispersion of collective modes because an expansion around a small four-momentum 
around the pole of K^^ gives the dispersion of the Goldstone boson. We suggest that it is the 
contribution of the Goldstone boson that leads the kernel to respect the gauge invariance 
condition. After making the analytical continuation i^i — )■ + ^5, the dispersion of the 
collective mode at T = is evaluated wtih — )■ and q — )■ 0. Note that q — )■ does not 
imply that the system is in the nonrelativistic limit. One necessary condition for that limit 
is hp ^ m. 



A. Nonrelativistic Limit 

In this limit, all contributions from the negative energy states vanish since A_(p) ^ 0. 
Therefore, we have 

p ^p+q^p ^ - (^p+q + ) 2ep+qep 

where -B(p, q) = ep_|_qep + Cp + p ■ q is defined in AppendixjF] and we have omitted the term 
proportional to v4(p, q) since limq_>.o v4(p, q) = 0. Due to the particle-hole symmetry, Qyi 
vanishes identically. Therefore Q'22 = Q22 so we only need to find a solution to Q22 = 0. 
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,2 



As q — )■ 0, ~ ^ — (/i — m) = 2m ~ ^ where /i~ plays the role of the conventional 
nonrelativistic chemical potential. 

At zero temperature, we keep the lowest order terms of u and q of Q22 and it becomes 



Q22(a;,q)=^ d^^ j dcosO- 



2 .Loo ^V£^p^ 3m^E^'^ 



(75) 



3 m 

where we have used = 2m(^p + /i^) and A^(0) is the density of states near the Fermi 
surface. Note that n ^ ep = in BCS theory so we found 

Q22(c^,q)^-^(c^ ^3^'' = "^^'^ ^^^^ 

where Cs = = is the sound speed of a BCS superfluid. Thus Q22(0, 0) = indeed 

and the expansion shows the dispersion u = Cgq of the massless Goldstone boson. 

B. Relativistic Limit 

Next we consider an ultra-relativistic BCS superfluid which is characterized hj kp ^ m 
and A <^ ~ ei? = ^y^^^-m^ ~ kp. Again the anti-fermion contribution can be safely 
ignored and due to the particle-hole symmetry we have Q12 = 0. Note that ^p+q — + q • 
Vf~ = + —■ Therefore we have 

1 , ,2 g^p^cos^e 
+00 pI — ^ ^ ^ 1 

2 / j^r^^a 



Q22{u;, q) ^ dpp J ^ dcose- 



47r2 Jo VEp3 3 el E-^ 

where = ep — has been used. The second term in the big bracket is at most of the 
leading order of Since m kp fi, it can be ignored. Therefore (522(0, 0) = and the 
expansion of (522(1^, q) leads to u; = -^g so the contribution from the massless Goldstone 
boson is clearly demonstrated. The dispersion also implies Cg = which is a well-known 
result for ultra-relativistic BCS theory. 

VII. COMPRESSIBILITY FROM RESPONSE FUNCTION 

The isothermal compressibility is given by k = n~'^{dn/dfi). Here the density suscepti- 
bility can be inferred from the response functions by [20|] 

|^ = -ir°°(a; = 0,q^O). (78) 
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At T = 0, the number equation (fT^ and the gap equation (fTSi) are 



n 



i = VC— —\ 



(79) 
(80) 



The density susceptibihty can be obtained from these two equations. We treat A as a 
function A(/i) of /i. Differentiating the gap equation with respect to /i one obtains 



(81) 



Differentiating the number equation with respect to /i and using the above result, one gets 

2 



dn 



(82) 



Next we check if ii''^°(0,q — )■ 0) can give the same density susceptibihty. When u = 0, 
Q23 — Q32 — Qi'2 — — 0- Therefore from Eq. (l36!) we have 



i^°°(o,q) = g;]!;(o,q 



00/ 



Q?3(0,q)Q!]i(0,q) 



Qii(0,q) 

In the hmit q — )■ 0, one gets -B(p, 0) = 2 and A(p, 0) = 0. Then at T = one has 



(83) 



Q-(0,q^O) = -2A^$:(^ + ^ 



Q?3(0, q ^ 0) = g°i(0, q ^ 0) = -2A 5^ ( 



Qii(0,q^ 0) 



gii(0,q^O) 



+3 



E-^ E+^ 



(84) 



After comparing this with Eq. (152]) . one finds that Eq. ( 175]) is indeed satisfied. This consis- 
tency implies that the approximation of the fermion propagator is compatible with that of 
the response functions. Thus we emphasize that the collective-mode contribution is impor- 
tant in maintaining the integrity of the CFOP formalism. 
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VIII. MEISSNER EFFECT AND SUPERFLUID DENSITY 



The Meissner effect can be demonstrated by examining the behavior of the response 
kernel K^^{0, q) as q ^ 0. In the previous section, we have learned that K^^ — Kq + 5K'^\ 



where Kq = and 



SK"-^ = Q11Q32Q23 + Q22Q3iQi3 - 2Q12Q31Q23 ^gg^ 

Q11Q22 — Q12Q21 

denoting the contribution from collective modes. However, in our model we found that col- 
lective mode effects do not contribute to the transverse components of the response functions 
in this limit and should not affect the Meissner effect. This is verified as the following. A ten- 
sor P*-' can always be decomposed into the longitudinal and the transverse parts Pl and Pt, 

where Pl = q--P-q and Pt — (X^i P^^—Pl)/'^- Assuming that q is parallel to the z-axis, in the 
limit q — )■ among all components of the response functions only Qg^^ and do not vanish 
to the first order of q. From this we conclude that limq^o Q3i • Q3j = hmq^o Q ■ Q3iQ3j • q. 
This means that the transverse component of the tensor Q3iQ3, vanishes in the limit q — > 0. 

Therefore the transverse part of K receives no contribution from the collective modes so we 
only need to consider K^^ — Q33 in the study of the Meissner effect. Defining p = p/cp, the 
longitudinal and transverse parts of the response functions are given by 

OL 



lim K,,, , „, ^ , , 
q_,0 J (27r) 



p p 



Now we focus on the nonrelativistic BCS hmit where ^+ ~ E^ ~ 2m (since // ~ m), 

Cp ~ m and fJ^~,A ^ m. Here a suitable rcgularization is needed to give physical results. 
This is done by subtracting the vacuum contribution from i.e., K^^[Q) K'^^{Q) — 
K'^ {Q)\T=A=o,i^=m- Finally we have 



2^n™ + 0(£j), (89) 



m m 

where is the fermion number in the nonrelativistic limit. Therefore, from J* = K'^^Aj 
we get the well-known London equation 



dE- 



m \ Sn^mn^^ .In dE' 
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--A(q)n„ (90) 

m 



where 

1 r°°.__4. df{E-) 



n, = n 



37r2mnNR L ^ dE, 



dPP'i--^) (91) 



p 



is the nonrelativistic superfluid density. Again the extra factor 2 in the London equation 
comes from the fact that we introduce the pseudo-spin a =t, J, so there are two times more 
degrees of freedom. 



IX. NAMBU'S INTEGRAL EQUATION FOR RELATIVISTIC BCS 

SUPERFLUIDS 

Before closing our discussions on the CFOP theory of relativistic BCS superfluids, we 
present a generahzation of Nambu's integral-equation approach. For nonrelativistic BCS 
superfluids, the spontaneously broken U{1) symmetry can be restored in the linear response 
theory by Nambu's approach ^]. In conventional BCS theory, the self-energy of the fermion 
is approximated by an integral equation which consists of a ladder approximation for the 
electron-phonon interaction. Nambu proposed that the EM vertex in the linear response 
theory should be corrected in the same way as the self energy. Hence the EM vertex function 
also follows the integral equation 

J2 a,G^''{K + Q)t''^\K + Q, K)G''^{K)a„ (92) 

K 

where the superscript "NR" denotes the corresponding nonrelativistic quantities. Explicitly, 
the solution to this equation should be an EM vertex that respects the GWIs. For the 
relativistic BCS model, following similar arguments we will derive the corresponding integral 
equation for the EM vertex. 

In Section IIIIl we found that the bare polarization function Kq"^ = (^33 does not satisfy 
the gauge-invariance condition q^K^^ = 0. This is because the collective modes which 
correspond to the fluctuations of the order parameter also contribute to response functions. 
The EM vertex without collective-mode effects in Nambu space is given by 

r = asY=[^;; 4.). (93) 
Then the bare polarization function (see Eq. fl33|) ) can be written as 

Ki;^{Q)=TTj2{rGiP + Q,f^)rGiP,fi)). (94) 
p 

The violation of the conservation law can be traced back to the use of the full fermion 
propagator and the bare EM vertex simultaneously [l|,0] so the approximations of the fermion 
propagator and EM vertex are not treated on equal footing. In quantum electrodynamics, 
the gauge invariance, or equivalently the Ward identity, can be maintained order by order. 
However, the generalized Ward identity is not respected even at the tree level for both 
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relativistic and nonrelativistic BCS models if the approximations for the vertex and the self 
energy are different. 

For the bare EM vertex in the relativistic BCS theory, one has 

q,r{Q) = ^^3 = crsGo'iP + Q,fi)- G,\P,fi)a3. (95) 

Note that the bare propagator appears in the identity. Thus a gauge invariant EM vertex 
must satisfy 

q^T^{Q) = a^G-\P + Q,fi)-G-\P,fi)as. (96) 

If we define the correction of the EM vertex as T^{Q) — 7^(Q) = Sr'^{Q), then from the 
two equations above and Eqs. fl59]) . one can get the GWI which is associated with the self 
energy: 

q^6t^iQ) = Sag - agS = 2tAa2tl5- (97) 



One possible method to respect the GWI flUBl) or flU7|) . as pointed out by Nambu |2|, is to 
treat the full EM vertex F'^ in the same way as how the self energy is approximated. That 
is, the full EM vertex of relativistic BCS superfiuids should satisfy the integral equation 



f^^Q) = YiQ) - 2(7 5^ 0-3 [g(P, fi)t^iQ)G{P + Q,^i) + G(P, -fi)t^{Q)G{P + Q, -/i) 



0-3 



riQ) -^aY^Yl ^^^)^'iQ)G{P + Q, cjii)a,. (98) 



P (7=±1 



To prove the gauge invariance of the above vertex V^{Q), we substitute Eq. fl^ into the 
GWI (1961) . After rearranging both sides, we only need to prove that 

2zAa2^75 = -2^7 5^ 5^ ^3G'(P, a^i)q,t^^{Q)G{P + Q, a^l)a^. (99) 

P <T=±l 

Details of the proof of this equation is in Appendix|Dl Therefore, the vertex given by the 
integral equation fl98l) respects the GWI, or in other words, the theory is gauge invariant. 
Moreover, from the proof in Appendix|D] we conclude that any truncation of the integral 
equation can not produce a gauge invariant vertex since terms of different orders of g in 
Eq. flD4p cancel each other. 

Interestingly, as pointed out also by Nambu j2|, for the nonrelativistic BCS theory, the 
integral equation is not only consistent with the generalized Ward identity associated with 
the EM vertex but also consist with the GWIs associated with three other interaction vertices 
(as shown in Eq. (4.4) of ref. ^]). Moreover, the integral equation of the EM vertex is a 
vector equation while the GWI is a scalar equation, they have different degrees of freedom so 
there is no strict one-to-one correspondence between the solutions to the integral equation 
and the EM vertex respecting the GWI. 

For the relativistic BCS theory, those conclusions should remain the same. Hence the 
integral equation fl98|) is not equivalent to the GWI associated with the EM interaction. 
What is equivalent to the GWI is the contracted integral equation given by 

q,t^{Q) = q,nQ)-2gJ2 E ^3G(P, (T/i)g/^(g)G(P + Q, a/i)a3. (100) 

P <T=±1 
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Since Eq. ( p8|l satisfies the GWI, we can also derive the GWI from Eq. fllOOp . Thus any vertex 
obeying the GWI (p6|) must satisfy Eq. fllOOp but not necessarily Eq. (!98|) . Substituting 
Eq. (!96|) into Eq. fllOOp . the right-hand side becomes 

P cr=±l 



P a=±l 



g^f'^(Q) (left-hand side). 



(101) 



Therefore it is the contracted integral equation f llOOp that is equivalent to the GWI but 
not the integral equation fl^S]) . A Comparison with the vertex determined by the latter 
shows that the vertex T''^ determined by Eq. ( llOOp may differ by a gauge transformation 
r'^ = F'^ + x^, where satisfies the Lorentz equation g^x'^ = 0- Such gauge transformations 
correspond to, for example, = 9^/, where / is a matrix of harmonic functions in Nambu 
space. 

In Section IIVI we derived the generalized vertex function S in the generalized gauge 
space. In Nambu space we would like to investigate the vertex T'^ and show that the GWI 
( |49|) will reduce to Eq. ( !96l) . Eqs. ( !34l) can be written as 



12 



Qi3 Q 
Q23 Q22 



Qn Q12 
Q21 Q22 



A,, 



Qi3 

Q23 Q21 



Qn Q12 
Q21 Q22 



A„ 



(102) 



where n^_2 satisfy Qfjli = and q^Il2 



-2iA by noting Eq. (l371) . Therefore the gauge 



invariant vertex F'^ is given by 

which obviously obeys the GWI 

q^T'^ = <^as - 2A(T275 = crsG-\P + Q,I2)- G'\P, 



(103) 
(104) 



As we have discussed previously, since f''^ satisfies the GWI, it should obey Eq. fllOOD . Hence 
it can differ from F^ given by Eq. ( |98l) by a matrix function xf at most. Moreover, the EM 
response kernel is now expressed as 



K^^{Q) = Tr ^ {V'^G{P + Q, /i)^ G(P, /x)) , 



(105) 



where Eq. (l36|) has been used. Compare to the bare response kernel (jUj), the GWI (I104p 
leads to the gauge invariance condition of the full response kernel qfj,K^'^{Q) = 0. 

Although a solution to the integral equation gives a gauge invariant vertex, it is a great 
challenge to find a solution. We emphasize that the integral equation (!98ll for the relativistic 
BCS theory should be implemented when one seeks gauge-invariant response functions. 
Previous work based on RPA approximations [sl implemented an iterative method without 
explicitly showing the complete integral equation for the relativistic BCS model. Further 
investigations of Nambu's integral-equation approach are needed for a better comparison 
between the results from the two approaches. 
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X. CONCLUSION 

The CFOP approach to the hnear response functions of relativistic BCS superfluids re- 
stores the gauge invariance of the response functions of an external EM field. The man- 
ageable computability of this theory allows one to explore several interesting phenomena 
including collective modes, compressibility, and Meissner effect. Importantly, this approach 
leads to a consistent expression for the compressibility. When the pairing interaction is 
tunable, a BCS superfluid may exhibit a BCS - Bose-Einstein condensation (BEC) - rela- 
tivistic BEC crossover [l2|. Interesting issues may be raised in the linear response theory 
of a relativistic BCS superfluid in the crossover. Our CFOP approach could provide some 
useful tools for investigating those issues. 

CCC acknowledges the support of the U. S. Department of Energy through the LANL/LDRD 
Program. 



Appendix A: Spinor conventions 

Here we use the Weyl or chiral representation of the 7— matrices. 



The metric is chosen as rj^i, = diag(l, — 1, — 1, — 1). The charge conjugation matrix C is 
defined as C = i-f'^-f^, which satisfies = -/, = = -C and [75,6"] = 0. The 
7— matrices satisfy 7^^* = and C'j^'^C = 7^". The energy projectors satisfies 

7°A±(P) = A±(-P)7°, 
i75A±(p) = A^(-p)i75, 

7V5A±(p) = At(p)7°^75, (A2) 



and 



(T+A_(p) + (T_A_(p) + a+A+(p) + a+A_(p) = 1. (A3) 



Appendix B: The Fermion Propagator in Nambu Space 

From Eq. ffTSj) . the Fermion propagator in Nambu space is 

G{P, /i) = [7°(ia;„ + f^cTs - 7°7 ■ P - "^7° + Ai7°75cri 



-1 

-1 



= [{iUnf - (A;^ + + + /x^) + 2^7°(m + 7 ■ p)a3] " 

X {iun - jJ^CTs + 7°7 ■ p + m7° - Az7°750-i)7°. (Bl) 

This can be written in a more compact form as 

G{P,fi)= (B2) 
[{iUnY - - - _ 2^/(m + 7 . p)a3] [iUn - /icrs + 7°("^ + 7 " P) - Ai7°75(Ti]7° 
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Note that -y^{m + 7 ■ p) = ep(A_,.(p) — A_(p)) and (13 = (T_(_ — o"_. The first part of the 
numerator can be evaluated as 



-el-A^-fi^- 2/i7°(m + 7 ■ p)a3 
= {{tu^r - E;') (a+A+(p) + a_A_(p)) + {{tcorf - E^') ((T_A+(p) + a+A^p)) .(B3) 

The second part of the numerator can be evaluated as 



lUJr, 



yWCTs + 7°("^ + 7 ■ P) - Ai7°75cri 



{iUn + ^p)(T+A+(p) + {iUn + ^+)ct_A+(p) + - ^+)a+A_(p) + (zw„ - ^p)a_A_(p). 



Note that 6 



±2 



.,±2 



)i?p and Mp^ + t;p^ = 1. One has 



p 

iUn - /ida + 7°(m + 7 ■ p) - Ai7°75ai 
= [Mp^(iu;„ + Ep") + v~'^{iuJn - ^p )]a+A+(p) 
+ [u+\iun + E;) + v;\tu^ - E+)]a_A+(p) 

+ b 

+ [ 



|Mp^(«CJn 



Up [^(^n 



p 

E+)+<2(*a;„ + E+)]a+A_(p) 



iw„ + Ep) a_A_(p). 



The denominator becomes 



(( 



4 - A' - /^') - V4 = ii^^nY - E^^mco^r - e; 



+2^ 



(B4) 



(B5) 



After substituting Eqs. ( ]B3[) . ( 1B4[) and ( IBSj) into the expression of G{P, fi) and using ct+cti 
(T+ and a_(Ti = cr_, we have 



E- 



+ 



^'p^ ^A / A , ^ "p^ 



+ 

+ 



u2 



+ 



^A / ^ , / ^p^ 



,-2 



za;„ - E+ iUn + E. 



-)A+(p) + ( 



+ 



iUn + E^ iUn - E. 



-)A-(P) 
-)A-(P) 



7°a+ 



7°a_ 



A+(p)A 



+ 



A_(p)A 



A+(p)A 



E+^ 



+ 



A_(p)A 



-E- 



ijncr-. (B6) 



Appendix C: (5-function operator and general properties of functions of e — Ep 

We would like to evaluate an arbitrary function with the argument e — Ep, where Ep = 
7°(7-p + m) — fids — A'y^ij^cTi. Following the derivation of Eq. ( lB4p (/x — i- — =^ E^ ^ E~ 
and 7°(7 ■ p + m) — )■ —7^(7 • p + m) ^ A_,_(p) ^ A_(p)), we have 

t-Ep 

= (6 + C)^+A_(p) + (e + )^-A_(p) + (e - )^+A+(p) + (e - C)a_A+(p) 

+ A7°Z75Cri 

= [u;\e + E;) + - E+)]a+A_(p) + [^-^(e + Ep ) + v^\e - E^)\a^K^{p) 

+ [u-^\e-E^)+v^\e + E^f^a^K+{Y^)+[ul\e-E;)+v;\e + E;f^a 
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+ u;v;{{e + E;) - (e - E;)) (a+A-(p) + ^-A+(p))7°*75^xi 
+ Mpt;- ((e + ) - (e - )) (a+A+(p) + a_A_(p))7°275ai. 

Explicitly, the four components of e — E^ in Nambu space are 

e - ^p = Up (e - E-) +v-{e + E^) + n+(e + E^) + {;+(e - E^), 

where 



Wp^A+(p) -npi;pA+(p)70i75 
-MpWpA_(p)7°i75 t;p2A_(p) 

^p^A+(p) Mpi;pA+(p)7°i75 
Wp^pA-(p)7°^75 np2A_(p) 



^A_(p) M+t;+A_(p)7S75 



<^^A+(p)7S75 t;+'A+(p; 

t;+2A_(p) -M+t;+A_(p)7S75 
-n+i;+A+(p)7°i75 



The four components of E-p in Nambu space explicitly are E-r 



\u. 



-2 



[v--u+^)e;k^ , ,.p .p 

-2(E+M+^;+A+ + E;Mpi;p A_)7S75 



t;p2)Ep A+ -2(i5;>+t;+A_ + E^u^v^ k^)^H^, 
'K" - <')^P^A+ + (^p2 - ..p2)i?p A_ 



(CI) 

(C2) 

(C3) 
(C4) 
(C5) 
(C6) 



.(C7) 



We have omitted the argument p of K±. After comparing those expressions, one finds that 

{E- + 4)A+(p) _ {E- - ^p)A+(p) 



u' 



2E- 



2E- 



(ig+-^p)A_(p) 
2E+ 



(i?+ + Ep)A_(p) 



2E+ 



It can be verified that the operators u and v also satisfy the following properties 



u 



±2 



(C8) 
(C9) 

(CIO) 



UpV^ = v^u^ = u^u^ = up^ = v^u^ = v^v^ = 0. (Cll) 
After multiplying e — £'p by itself, one gets 

(e - EpT = S - ^P )" + ^P + ^P )"■ + ^P + ^P )" + ^P - ^P (C12) 

where n is a positive integer. One can find that Eq. flC12p holds for n = 0, too. Interestingly, 
from Eq. flB6p one can see that the case n = — 1 is also valid by inspecting the expression of 
G(P,/i)7°: 

(e - EpY' =Up{e- E^r' + Vp{e + E^y' + n+(6 + E^y' + v;{e - E;)-\ (C13) 
Following the same argument, 

- ^p)-" = S (e - Ep)-- + Vp{e + Ep)- + n+(e + Ep+)-" + t)+(6 - Ep+)-". (CM) 
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That means that Eq. ( 1C12I) holds for any integer n. For any function F{e — Ep) , we have 
the expansion 

Fie - 4) = WpF(e - E^) + v^F{e + E^) + u^Fie + + {;+F(e - E;). (C15) 

Thus, 

e-«(e-i?p)* ^ ^-e-*(-^p )* + {)-e-*(^+^p + n+e-'(^+^?)* + t)+e-'(^-^?)*. (C16) 
By Fourier transform we get 

5{e - Ep) = u-J{e -E^) + v^6ie + E^) + ^+5(e + + {)+5(e - E^). (C17) 



Appendix D: Integral equation of EM vertex and GWI 

Before proving that a vertex determined by the integral equation (l98i) must obey the 
Ward identity (jH]), or equivalently, verifying Eq. fl^ . we prove the following identity. 

'^9J2Y1 ^z") - ^(^' ^'")^3) = -22Aa2^75 = S^s - ^sS- (Dl) 

P o-=±l 

The left-hand side of (IDip is 

Ep(F(P,/x) + F(P,-/i)) 

-2.pl 2i?- + 2E+ nS 

-2zA(T2^75, (D2) 



4^? 
4zA5( 



where the gap equation (fT5|) has been used. From E = Gq ^(P, /i) — G ^(P, /x) one concludes 
that 

G{P, fi)Go\P, fi) = l + G(P,/x)S, Go\P, fi)G{P, fi) = l + SG(P,/i). (D3) 



Now we turn to the proof of Eq. fl99l) by considering the right-hand side. By substituting 
Eq. fl98l) into the right-hand side and repeating the process, we get an iterative equation 

RHS of Eq. (EH]) 
= -2g J2 asG,{P)q,riQ)GAP + Q)^z 

Pa 



+ {2gf ^3a(Pi)a3G'.(P2)g^7^(Q)G.(P2 + Q)a3a(Pi + Q)a; + ■ ■ ■ 

OO 

Y,i-2gY Yl ^3a(Pi)---a3a(P0g^7'^(Q)G.(P. + Q)a3---a(Pi + Q)(T3 



= 1 Pl-PilT 



OO 



1=1 Pi-Piak=l k=l 
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where we have defined Gcr{P) = G{P,afi) to shorten the expression. After inserting the 
Ward identity (!95|) for the bare EM vertex and using Eqs. (]D3| ). we get 

RHS of Eq.(l99]) 

oo 

= J2i-29y <y^G,{Pi)---a;G^{Pi)a:iG^^{P, + Q)G,{P, + Q)a:i---G,{Pi + Q)a:, 

i=l Pl-Pi<T 



OO 



-J2{-2gy (rsG^iPi)---asG,iP,)G,,\P,)asG,{P, + Q)a,---G,iPi + Q)as 

i=l Pi-Picr 

= -2gY{^3G^iP) - G.iP)as) 

Pa 

oo 

+ Yi-'29y Y ^3a(^i)---^3G.(p.)a(p.-i + g)a3---a(Pi + g)a3 

1=2 Pl-Pi(T 

oo 

+ J](-2(^r Y (T,G^{P^)■■■a■,G,iP,)a,±Ga{P^ + Q)(T^■■■G,{P, + Q)a, 

i=l Pi-Pi<T 
oo 

-Y{-2gY J2 cr3G^{Pi)---cxsG^{P,-i)G^{P, + Q)crs---G^iPi + Q)a^ 

1=2 P\---Pi(T 
oo 

-^(-2(7)* (r^G,{P,)---a,G„{Pi)ta,G,{P, + Q)a^---G„{P, + Q)a^. (D5) 

i=l P\---Pi(T 

By changing the dummy index i — )■ i + 1 for the second and fourth summations, we get 

RHS of Eq.(l99]) 

= 2iAcr2275 

oo j i 

- Y^-'^SY Y n [^^G.{Pk)\2g Y [cr3a(P.+i) - G.(P.+i)a3] H [G.{P^+l-k)cT,] 

1=1 Pi-Piuk=l Pi+i k=l 

oo i i 

+ Y^-^aY Y n hG^iP'^)] - Eas] J] [^.(^.+1-^)^3] 

i=l Pi-Piak=l k=l 

= 2iAa2275 = left-hand-side of Eq.([99]), (D6) 

where Eq. (]Dip has been used. Therefore we have proved that any vertex that satisfies the 
integral equation must also satisfy the Ward identity and hence must be gauge invariant. 



Appendix E: the odevities and symmetries of the response functions 

Since the energy projectors satisfy 75A±(p) = A=p(— p)75, we have (See Eqs. (fT2i) and 

m) 

75G'(P,-/i) = G(-P,/i)75, 

75P(P,-/x) = P(-P,/i)75. (El) 
The energy projectors also satisfy 7oA±(p) = A±(— p)7o so we can conclude that 

7°G(P,/.) = G(P,/.)7°, 
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7°F(P,/x) = -F(P,/x)7°, (E2) 
where P = {ioOn, —p)- The propagator in Nambu space flB6P can be written as 

G{P, fi) = G{P, fi)a+ + G{P, -fi)a- + F{P, fi)a+ + F{P, -fi)a.. (E3) 
Using Eqs. (lEip . a2cr± = —o^a^ and 0"2(t± = o^a^ one can show that 

(T275G(P, ii) = (r-iG{-P, ij)a^a2i^. (E4) 
Similarly, using Eqs. (lE2p . cr3cr± = — ct±ct3, and a^a± = (T±(T3 we have 

a37°G(P,/i)=G'(P,/i)a37°. (E5) 
Thus we can analyze the odevity of the response functions of the four-momentum. Explicitly, 

Q,,(zfii,q) = ^Tr[(o-275)'S.G'(P + g,/i)S,G(P,/i)] 
p 

^ ^_^)(2^-+25i.+5-+5^.) J2 Tr [±M-P - Q, /i)S,G(-P, /i)] 

p 

= (-l)(^"+^^^)^Tr[S,G(-P-Q,/i)S,G(-P,/i)]. (E6) 
p 

Changing variables by P — )■ — P, one gets 

Q,jitni,q) = (-l)(^"+'5'^)^Tr[S,G(P-Q,/i)S,G(P,/i)]. (E7) 

p 

Therefore, 

Q,,(zn,,q) ^Q,,(g) = (-l)(^"+^'^)g,^.(-Q) ^ (-l)(^"+^^^)g^^.(-,f]^,_q). (E8) 

Using the relation (]E7p . one can show that 

g,,(zfi,, q) = (-l)(^"+^'^)Q,,(zfi,, q). (E9) 

Next we analyze the odevity of the response functions about the spatial components of the 
momentum. Using Eq. (lESp . for i,j ^ 3, one has 



Q,,(zfi,,q) = ^Tr[(a37°)'S.G'(P + Q,/i)S,G'(P,/i) 
p 

= 5^ Tr [t,G{P + Q, fi)t,G{P, fi)] 
p 

= 5^Tr[S,G(P + g,/i)E,G(P,/i)] 
p 

= Qi,itni,-<i), (ElO) 

where in the third line we have changed variables by P — ?■ P. For i = 1, 2; j = 3, with the 
help of Eq. (]E9P , we only need to consider the case with i < j. Thus 

QU^^i, q) = i-iy^'^^'Q^si^ni, -q). (Ell) 
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For i = j = 3, we get 

Qt;^,{tn,, q) = i-ir'^''"Qt;^,itni, -q). (E12) 

From the odevities of the response functions about the four- momentum (see Eq. (]E8|) ) and 
the spatial momentum (see Eqs. (lElOp . ( lElip . and f lE12p ). we derive the odevity of the 



response functions about the boson Matubara frequency as follows. For i = j = 1, 2 we have 

Qii{iQi, q) = Qii{-iQi, q). (E13) 

For i = J = 3 we have 

Q'ssm, q) = ( ^BT'^^ ^\ -/^ ^ = " = ' " ^ = = \ (E14) 
"■^ ^ ' \ —Q^J,\—^^lh^) II = 0, = z or /i = z, z/ = ^ ' 



For i = 1, j = 2 we have 



Qi2mu^) = -Qi2{-i^uc^)- (E15) 



For z = 1, j = 3 we have 



For i = 2, j = 3 we have 



Appendix F: Expressions of the Coherence Coefficients 

For convenience, we introduce = (ep, p) and = = (ep, — p) so the energy 
projectors can be rewritten as 

, ,, ^ + 171 r. nt + m 

2ep ^ ^ 2ep ' ^ > 



which also satisfy 



A+(P)75 = 75^^ 7 = 757 = 75A-(-p), 

ZCp ^Cp 

A-(P)75 = 75^7° = 757°^ = 75A+(-p) (F2) 

We define A(p, q) = Cp+qCp — Cp — P ■ q and -B(p, q) = Cp+qep + Cp + p ■ q, which satisfy 
A{'Pi q) = ^(— P, — q) and -B(p, q) = i?(— p, — q). By using the above relations and the the 
identity 

Tr [{^' + m)7^(|^ + m)7^] = 4 [A;'^^ + A;"^F - ^^^(A; ■ A;' - m^)] , (F3) 



we can show that 

(u u )i^ = {v V )ii = ^(l zr- F- )- 

{U U )22 = {V V )22 = :^(1 )- — . 

^ -^P+q-^P ^P+q^P 

If// = i/ = Oor/i = i, z/ = j, the 33-component is given by 

( 1 + Meis) if ^ ^ ^ 

(^1 + gp+.gp+A^ A (p+q)^p^- + (p+qyp'+^^^A(p,q) ^ ^ ^nd = 

V ^p+q^P / ^P+1^P ^ 

lf/x = 0, u — i or II — i, 1/ — 0, the 33-component is given by 

^ ^-^P -^P+q^ ^P+qfp 



{U M )?3 + (v W ) 







A(e- + e;+q)^(p,q) 

^p+q tp+qtp 



^ ^p+q ^p ^ ^p+qf p 



(ii-ii-)°3 - (^;-^;-)°3 = ^a(-^ - — ) 



1 \i?(p,q) 



{u u )\^ + {v V Y 
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p+q -'^p ^p+q^p 
^p ^p+q £p+q^p 



(i^ = u )ii = -(1+ r^- )- — • 

-'^P+q^ ^ ^p+q^p 

(li ^; )22 = (^^ )22 = :^(1+ ) • 

^ ^ -^P+q-^P ^ ep+q*^p 

If/x = i/ = Oor/x = i,i/ = j, the 33-component is given by 

^p+q^p / ^p+q^p ^ 

_ gp+qgp+^' \ (p+q)'p^' + (p+qyp'+y^A(p,q) jf ^ ^ 

•Bp+qEp y ep+qep ^ 

lf/x = 0, 1/ — i or II — i, 1/ — 0, the 33-component is given by 

^-^P+q ^P ^ Cp+qCp 
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(u V )i2 = -{v u )i2 = + — • (F16) 

(u'v-)UM,-^%:^^, (FIT) 

^-^P+q -^P ^ Cp+q^^P 

(^^-.-)°3 - {v-u~)% = zA(-^ + ^)^^, (F19) 

^-^p+q % ^ ^P+q^P 

-'^p -'^p+q ^p+q^p 

{uV)n^{v^v+)u^ -il ^+ p+ ) — ■ (F21) 

^ ^ -^P+q^P ^ ^P+q^P 

= = i (l - (F22) 

^ ^ -^P+q^P ^ ^P+q^P 

If = 1/ = or = i, u = j, the 33-component is given by 
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f 1 + ^hrlf ) ^^^^^^^^'^^"-'^^^^ if ^ i and ^ ^ j ^ ^ 



^p+q^P 

11^ — u — Qor^ — — the 33-component is given by 



— ^p+q I ^ ^p+qP + ^p(P + q) (-{?oA\ 

[U U - [V V - - I — p h ^ I . [t 24 j 

^^P+q ^P^ ^P+q^P 



2 



f ^P _ C+q \ -B(p, q) 



p 



p+q ^p+q^p 



+ («+«+);3 = '"g,t5"'''f^'" ■ (F26) 

-C'p -^p+q f p+qf p 
^-^p+q ^p^ ep+qCp 

(z.+z.^)°3 - (.^.^)^3 = -^a(-^ - ^) f^^, (F28) 

V-t/p+q -i^jT^ ^^P+q^P 

(.^.^)-3 + (.^.^);3 - Zf*"' + ^ P'^^^r (F29) 

J^p ^p+q Cp+qtp 
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P+q P "^p+q'^p 



Now we evaluate the "mixed" terms. 



/ - +^ / - +^ ^P+qC + ^' A -4(p,q) 

^ ^ ^P+q^P ^ ^p+q^p 

/ - +N f - +^ Vi C+q^P - -4(p, q) 

(ti t,+)22 = ^;+)22 = - M j — 

^ ^ -^P+q-^P ^ ^P+q^P 

li II — u — or II — i, u — j, the 33-component is given by 
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+ _i_A2^ 



-fl + ^'^:^^'^:f ') (p+q)'p-+(p+q)-p'-^'-^(p.q) if = i and = j. 



(F30) 



. + +X / + +^ Vi , ^P+q^p +A ^ ^ P,q .-pon 

^ ^ -^P+q-^P ^ ^P+q^P 

If// = i/ = Oor li = i, = j, the 33-component is given by 

fl_%^)MP:H) if;. = ^ = 

f 1 - ^hrtf ) (P+q)^P^+(P+q)^P'+^^^^(P-q) if ^ i and ^ ^ J ^ ^ 

V K+'i^i J 'p+1'p 

lf/x = 0, 1/ — i OT n — i, u — 0, the 33-component is given by 

{u+v+)i2 = -(^^ V)i2 = I T^q: + ) — ■ (F34) 

^ ^-^P -^P+q^ ^P+q^P 

(u-v-)1, + (v*u-yi, = ^%±$^^, (F36) 

-^p+q ^p+q^p 



^-C/p+q -C'p^ Cp+qCp 



(K+.+)^3 - (.+«+)^3 = -^A(-^ + ^) (F37) 

. + , r„+„ + V _ ^^(^P - ^P+q) (P + q)'^P + P'^P+q /poo^ 

^p ^p+q Cp+qtp 



(F39) 
(F40) 



(F41) 
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li IJ, — 0, u — i or ijl — i, u — 0, the 33-component is given by 

{u-u^rss - {v-v^)^l = (1^ + ^P + q)^^p-P^^P+^ (F42) 



^P+q^ tp+qtp 



-^p+q-^p Cp+qCp 



(F44) 

);3-(.-.^);3--a(^-J^) 'p + ''>'-'-p'-^- . (F45) 



(^^-^.+)°3 - (.-.+)^3 = iA(-^ + ^) (F46) 



^P+q^P ^P+q^P 



( - +\ ( - +\ , ^P+q^P + A'^ A -^lp.q) , . 



p+q-^p ' "-p+q"-p 



(F49) 



(F50) 



^ ^ -^P+q-^P ^ ^P+q^P 

If// = i/ = 0or// = i,i/ = _7, the 33-component is given by 

fi- y^^+f )^ if^ = ^ = 

_fl _ ^P+'^^'^-f ) (P+q)'p-+(P+qyp'-^'-i^(P,q) if ^ = i and = J 

lf/x = 0, 1/ — i or fj, = i, u = 0, the 33-component is given by 

(^^-.-^)- - (.-^^+)- = - 1^ ) (P + - P^^P+^ (F51) 

^-C'p+q -C/p^ Cp+qCp 

(.-.-^),. = = + 1^) (F52) 

^ ^-^P -^p+q^ ^p+qep 

^p+(i ^p+q^p 
^^p+q ^P^ ep+qCp 
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{u-v+)l, - {v-u+)l, = tA( - J-)^l(Pl^, (F55) 
iu-v^)^. + iv-^t. - -^^§^^^±^^^^^. (F56) 

P+q P Cp+qCp 

r + ( + -\ Vi ^P+q^P + -4(p, q) 

(li+ii )ii = [v^v )ii = I 1 p+ p_ ) • (F57) 

^ ^ -^P+q-^P ^ ^p+q^p 

[u+u )22 = (^;^^; )22 = 7^ 1 1 p+ p_ ) - — — • (F58) 

^ ^ -^P+q-^P ^ ^P+q^P 

If/x = i/ = Oor// = i,i/ = j, the 33-component is given by 

+ if/^ = ^ = 

-fl + ^^^^' l (p+q)'p-+(p+qVp'-^'-^(p.q) if = i and = j 

lf/x = 0, u — i or n = i, u = 0, the 33-component is given by 

/ + -\u,u I + -xoi' ( ^p+q I \ (P + Q) ^p ~ P ^p+q l-c^cr\\ 

[yTu )^3 - (v^v )^3 = - + ) ■ (F60) 

^-^P+q -^P ^P+q^P 
/ + -^ ^ + V^P ^P+qA^(P'q) /T^fiiN 



(F59) 



(F62) 

M(p + q):!p^P:!p±q (pg3) 



13 ^ yi3 

^^P+q ^P ' Cp+qCp 

^ -'^p+q ' £p+q£p 

p+q P Cp+qtp 

(i^+t; )ii = [y^u )n = - 1 1 + „_ ) • (F66) 

^ ^ -^P+q-^P ' ^p+q^p 

= (V+Ii )22 = - I 1 + „_ j . (F67) 

^ ^ -^P+q-^P ^ ^P+q^P 

li II — ly — or /I — i, u — j, the 33-component is given by 



^33 
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-fl - (p+q)^P^-+(p+qyp^-^^^-B(p,q) if ^ = , and = J 



p+q 

If// = 0, 1/ — i or II — i, 1/ — 0, the 33-component is given by 

':+ 



(F68) 



(«+^-)^3^ - {v^u-r,^ = - - ^) (P + q)^^P P^^P+^ (F69) 

^-^P+q -^P ^ ^P+q^P 

(W+I; ),2 = -(v+M )^^ = - I +-^\ (F70) 

(„-.-)!3 + (.-«-)!3 = ^^r-^p-'fp-") , (F71) 

-^p+q-^p fp+qCp 

(«*«-);3 - = -a(^ + l) (P + q)Vp-p-.p., 

^^p+q ^p^ ep+qep 

(«+.-)^3 - (.+«-)^3 = -^A(-^ - ^) ^i^Hl^, (F73) 

^ -^p+q -^p ^ ^P+q^P 
p+q P Cp+qCp 
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